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: Thermodynante properties of all commonly. u used m. “Aittery gun. 


propellants. plus some experimental propellants are presentca in 


| tabular and graphical form. For the military propellants ? ‘these 
properties are given over a ‘loading density ranging from 0.05 ‘to 


0.4 en/ee. Compe: isons. are made between the propellanta on the | 
basis of "force" nd Laochoric flame temperature. The tables and - 


data presente t he: ein can be used for gun, dente computations... 
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: ‘Within the past two. years a general prograz for the computation —~ 
, of ballistic trajectories inside or a gan‘) has deen set up at these 


% ; latoratories on the high speed digital computer p CROVAC. A portion - 


of the input data’ re, muired for this program is thermodynamic prop- 


'. erties of the combastioa pret:cts such a3; ; flare vemperature, force, . 


Sern, aud specific keat ratio. 


Tables which have bead eed in ‘tae past es @ source of this 
“{efornition ‘have. been compiled either from obsolete *shermocynamic — 
. Gata, ‘os computed using various simplify ying assuaptions. In addi- 
ti un, these tables have given the therzodynamic properties for a single 


‘@annity of loading for the gun (usually 0.2 gfe) vhereas most guns 
Operate over a wide range of loading detsit! es. Since the thermo- 


dynar' ‘properties of the ‘propéllant gases var f with loading density, 
te 


2% would be useful to compute the properts @8 of all the common 


military propellants over the range of Aoadiag asxsities normally 
eneountered in gua firtngs. : 


“In ‘the following seport, the thernodynente properties for all the 
corecn mil{tary propellants and some experinental propellants are 
pice Pes thy military propellants (with the exception of M/ and 
Meo) shese propesttes are given over a loadiag density range from 
0.05 to 0.8 gafec. In at’ “tion the therzodynante properties of 
three other types af prow ianta are exaxined: 


(2) Nitrocellulose at verying percontages of nitration. 
(2) ww propellant contaalaated elth varyiag percentages of © 
watex and/or ethyl alcohol. 


(3) Double dbase propellant vith saviae ere;ortfons of altro- 
cellulose ond nitroglycerine. 


To carry out the extensive computations reguired for the pre- 
peretion of thie regort, a digital cocputer progras designed to 


i 


ae " compute the thermodynamic proverties of gun propellants vas used on 


ee ee ee ee This progran 


was initially formlated by Brinkley‘? ). 
on # and revised by Leger'3 ) end has been used by the authors for about 
“Oe ee ea : 1/2 year, The computer program as devisea by Leser has been revised 
e = a Me, nak ior aaee by one of the authors (Baer so that the flame temperature of the — 
ae Sat en Lak Re propellant can be computed together with the thermodynamic properties 
ae ‘A ol degree . gue ¢ as a av the flame temperature. This revised computer progran is the one 
eo NES 4 ae _Tow in current MRE e o ee bot 
bis | “IT COMPUTATIONAL RESULTS 

The data in this report have teen organized in the following 
ay Satie Ge a (LY Gables listing propellant composition and important thermo. — 
mS , ones es a a eynante ener eee a ak ere 


ov 


a Pe tay Tables of input thermodynamic data used in the program, — 
EE: -: (3) Figures presenting tn 


x : graphical ‘form some of the particular 
Ss) thermodynamic properties of military propellants. 
a ; (4) An appendix preacnting the detailed Computed data for all 
ats i propellants considered. 
ve vk 7 << oe For eese In locating & particular plece of detatled information 
as oe ; in the appendix, a table of contents appears at the beginning of that 
ne section, | 

78 In Table I, composition and thermodynamic data at the flame 
am ” temperature are elven’ for 16 atlitary propellants. In Table II, 

s . roe A _ data are given for T experinental opellants, In Toble IIT dota 

i BS : . are given for M1 prope Llant contaminated with water and/or ethyl 


alcohol ¢*4 for sitroce Liuloge propellants with varyieg degrees of 


altration. In Tudle IV data are xf on for experimental double base 
propellants, 


Tables fT through tf have heen arranged tn a Tortat slallar to 


-hat of the Picatlany tobles(* 


. ay abe 
eee: 


an 


vith &.. follovin; differences; 


So Ee 
ceils 
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(1). For éanh eee in Tables tr, It, ; ané ITT percentage 


* compositions: are given; the: composition under the column marked *s" 


(for specification) is the composition as given in the Picatinny 
tables. 


‘isa normalized composition used. in computing the thermodynamic 


‘properties of the propellant. The difference between the two. 


_ compositions is due to two factors: (1) percentages Listed in the 
Pleetinny tables added up to values greater than 100%, and so had 
_ to be normalised to a: percentage composition; (2). in sone propellants 


“ inorgante salts, such £3 potassium aulfate, petascius perchlorate, 


etc. hol been added. Since the Brinkley code does not have provision 


. for the tnelusion of these salts , Shey were ‘normalized out of the 
“ay composition. oars : 


It is believed ae with the exception of M7, the oalasion 
of those salts will not appreciably change the thermodynantc results. 


(2) Data in Tables I, II, end III are reported ot a density of 
“deadline of 0.2 gn/ce. 


loadiag of both 0.1 and 0.15. The difference in the densities of 
leading between Table IV ond the other three tables ia that the 


data in Table IV vere computed in the support of an experimental 


program for the measurement of ea gas temperatures in a 
' eloted bond. 


(3) The propellant Porce “PY tn al) four tables was computed at 


tite loading deasity shown. This waa done to give a direct compartson 


between the force values for the military propellants and the experi- 
ental propellanta. 


It te of interest to campare the force and flane tesperatures 


for these propellants. In Table V, aflitary end @: srimental pro- 


pellents are listed ta onder of decreasing feree; in Table VI ia 


order of decreasing adiabatic flac temperature, It will be noted 


that for some propellants there is a difference tatwren the force 


The composition under the column marked "Bp" (for Brinkley). - 


Data in Table IV are reported for densities of - 


Pear (ete aan 


rating and the flane tempersture rating. The special propellants of 
Tables Inf and IV were not included in this aio in ‘order to avoid 


ge confusion. - 


The thermodynamic data in table I are given for a standard losding hoe 
‘density of 0.2 gm/ce. To examine the variotion with loading density 


'". of the thermodynamic properties at the flase temperature, the thermo- 


dynamic properties for all the. military propellants (vith the exception - 

of NT and N26) were computed for five loadirg densities: 0.05, 0.1, 

' 0.2, 0.3, 0.4. This range was chosen to account for most of the gun 
. pressure range (1000 - 100,000 pst}. Table VIT gives these therno- a 
_ dynamic properties for the above loading densities, 


One of the results vhich can be derived from the Brinkley thermo- 
_dynanic computations is the covolume correction. This correctioa 
appears in the Nobel - Abel spetie of state for r propellant gas as 
follows: “ =. 

Ch a a) 

Since the voluse term is equation (1) is the reciprocal of the. 

\cading density, equstion (1) can be written as follows: 


| 3 -y) =F a) 
Expanding and rearrangirg (2) ve odtata; 
Pe FeyP. 3) 


Equation (3) can be represeated as a lincar plot of P/p ve P 
fa which the slope {s tie covolure ¥, and the intercept vith P « 0, 
the propellent force, F. Whea P/p ve P to plotted from the Brinkley 
computations, the resulting plot ts nat lfarer, the data pointe 
' falling on a smooth curve. Por this reason, the covolime carraction 
Ul) vary with the constant volute chanber ereasure P. To deter- 


nine thie varlation a cuble eqiation In pressure, P, was ritted by 
least squares to the date of Pye vs P, giving the following reitation: 


ae 


<f, a4 
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“The covolume is ‘ie first derivative of seuabion (3) thus: - 


2 = “8 + 2aP + ‘Sap (5). 


In equation (5) a, is the linear covolume tern. 
the value of the coefficients F, a,, 2a, and 4a 


283 are listed for the 
military propellants. Jn Table VII, for each military propeliant, 


_ the covolune correction is given for each density of loading. - 


In Table VIT is also given the intercept propellant force F, 
and a force computed from the relation: aa 
. 7 7 oo a 7 . : “ . 6 
F ann, . (6) 
The force F is tabulated for each propellant and for iach 
‘density of loading. 


The intercept force PF can be used to compute the pressure 
developed by burning the ee stiait fn an adiabatic closed boub -t 


@ loading density . “The pressure developed can be computed by 
solving the cubic equation: 


esr + pag +(ea, - 1p +¢F, = 0 (7) 

Since the value ‘for the coefficient ay ie quite small, Little 
error wild be made at low loading densities {f one drops the cubic 
term, in which case solution to the resulting quadratic equation 
can be obtained by norm) algebraic aethods. Coaplex algebraic 
solutions exist for the cubic equation, ose method being Msted 
in the Cheaieal Rubber Handbook. '? 

A nustcre of silitery guna, especially th-se of the hovitzer 
clase, can be approxiaated by a aisple gun adel. Ta thie sodel, it 
is assumed that the projectile does not move ti -ny appreciable 
exteat while the sropeNant {6 burning. After the propellant ta 
burned, tte projectile is accelerated by the exganding hot gaees. 
It fe of interest to toterlor ballistics to stov the relation 
Datween gas pressure and tenperature vhile the coxbustion gases 


ore undergoiag an (sentropte expansion oehlad the profectile. A 


In Table VIII 


Sortie «fy, 
mT v3 PES ee pa a? aa 


pressure-entropy dingram 43 one way in wt tch this process can ba. 


illustrated; therefor @; pressure- -entropy diagrams ave ‘been plotted 
from the results of the Rrinkley computation on the militury urOS 
pellants. These plots are chown in Fig. l-ik 


These ‘figures are a graphical representation of the relation- io 
ship. between vas density, temperature, Pressure, ay entropy ror the 
combustion gases of each ot ‘the military propellants, The dotted oe 
line labeled ™! abe aah the flame temperature of tne “propellants. 


These ftperes ean be u: red to predict the gas teriperature in a 
. gun, “provided that the pressures ore known. ‘For exumple ; suppone mM 
propellant at a density of londing of 0. 3 gmfec is ‘burned ina gun ey 
chanber andl burnout of the propellant occurs before. ‘the projectile 
. moves to ony appreciable: extent: . On Fig. l, the pressurs-cr ntropy 
dlagram for Mw propellant, the Latersect.on of the density of loni- 


ing Lue 0 « 5 with the iyochoric Plane ‘Lomperature Line tT, is 


the condition, in the chamber at. burnout of the propellant Under. 
these condi tons T, = 2440" XK, the press sure Pp = 60,700 psi... 
‘(this could: be a againsh thie experimental pressure) » aid the 
entropy 8 #. 2.122 cul/gm°K. If one ‘accuses that motion of. ‘the . 
"projectile causes inentro,ic expansion of ‘the gas; ‘then one can 
determine: the temperature of the gas-at any point along. ‘the ean 
tube provided the pressure at that point. is known For crample, 
assume that muzzle~pressure of the gun has been found to be. 5,000 

psi. Than dropping down the entropy line at 2.122 until the ‘5, 000. 
ae pressure Line le reached, one finds that the gas temperature is 
1540 K. are : 
. In closed chamber irings, the. mauke nohip between press sure and 
‘ohpetcurea Jf uenerally unknown because of heat 3088 from the 
chenber.s From these figures the temperature of the propellant gas 
cau be pr edicted, knowing the loiuling density and the experimental. 
pressure, ‘stneo the gas tempera, ure ile be the intersection of the. 
pr ensure dine with the was denaolty Line. 


4] 
rd 
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: Other ‘tyoas of aiagrans can be edectaact from the thernodycante 
. “date obtained from the. Brinkley couputations, but in “Meu of plotting . 
“more diagrams the thermodynamic: data is presented in the appendix vat 
for each propel lent. In addition to the therzodynanic anta, there | 
‘are presented the equilibri: combustion gas composition ata 
‘particular gas density and teuperature. : The ~ange cf gas temperatures: ae 
ned presented in these tables depends upon | the’ range oF ges ‘temperatures ie, ; 
, presented in the presstre-entropy figures. Fox high gas densities the | 
; ‘temperature Fenge is narrow. and for low gas densits ies, the temperature. 
: ; ‘range is wide. . Further details on these tables is presented in the: 
: appendix. . eS ; 
ee MISCELLANEOUS COMPUEATIONS, 


rt is of. interest to the: -nanutacturer ena user of military pro- : 2 
- pettants to deteraiie the ‘effect of water’ and/or ethyl alcohol on — ; 
o __ the therwodynantc properties (particularily the: flame tempiaeture, 
ee and: force) or the propellant. “The propellant chosen ‘for thie 
2" tavestigation vas M1 at & loading density ot 0.2 enfece the first. ; 
. zc portion of Table wit gives the propellant. composition and a cummary . = 
of the propellants therusdynamic ‘properties. . The format for this 
table is the came as that used in Tables I-and IL. The. propellant: " 
Composition column: headed by the letter "6" 45 the standard compoai<’. cae 
~ tion showing ‘the amount: ‘or alcohol and/ or water added, The composi: 
tion column headed by the letter "pt is the normalized aes ie 
a : composition used in the Brinkley computations... -In this table - 
amounts. of water. or ethyl alcohol added to ‘the propellant ‘were 
varied from 0 She ‘In one propellant a ip - 3 mix of ethyl _ 
alcohoi - ‘water’ Was added to determine the eft fect of a mixture. TA. 
plot of the flame temperature Vs percentage water or ethyl alcohol: 
- added is ‘shown ia Fig. 16. This figure can: be used to estimate the . 
effect on flame temperature of alcohol - ater mixtures , which can | ; 
be unite by eee ee the arop ia flaie temperature ne to water | - “ 


“: C ee 


oe ee a at PE! AIT CE EA 
et ee Ae Ve ae : : : 
ee ae ; , 

Ss a ne ean bé ‘adacd to the e drop due to. 5 alacbat “on the graph this estinate 
tae a te Ta "ean be made hy. measuring the temperature drop due to the addition of 

a é ee ee aie ia "water with ‘dividers and then apacing orf this. distance below the, $ 

4 ‘flare ‘temperature due to’ the alcohol edition. For tov (less. than ag) 

i. . . - percentages ‘of alcohol. and vater, this method will give values within’ 

e ; ee oe =. 3% of the computed value. - For higher percentayes, the error is: 

“ d , eae greater, for Anstance a 3p - 36 mixture gives on estimated vialue of 

ay PLT, the conputed ‘value being 2192° a Detailed. tabler for: hese ~ 

3 : _ Propellants are presented in the appendix, a eee ee 

 - an the Latter” portion of table paey are shown the. ebuibes teten eae 

Ss and thermdynamic properties of pure nitrocellulosé propellants where .. 

Bs ; a “nitrogen content varies from 12. oF to 13.0%. the Getatlea Somat 

“8 ee : , "tional results are + given in the append, a Bae ee ee at ated 
Be Recently, at this Yaboratory;, laine temperature abusias's were Se 
be ta Se nade on a series. -of double-base propellants - da which the nibtzon : 

is aoe se Se glyerine concentration was varied from Ob to 20p, Tm Table we 

a. Pee ie 28 Pe Se 


the: computed’ composition and thierinodyhamte properties, for ‘these. 


“experinental, double. base, Brops‘iants: ‘ave given for. two Joading ” me te Uae 

a Pe eer densities, 0.10 and 0.15 enfec. “As, ‘vould ve expected, ‘the flaise as ee 
ee i are ees temperature and ‘the’ propellant force increase as the aitros | 

re . PSTN ghd wlycerine content Aneeases Détatled conputations ‘on _theee - 

i te aoe eS ee ., Propellonto are elven ‘in the: appe Shain, 


_:- eee ae oe Th al’ of ‘the compittations, the combustion “products have 
1 RI SN Ser mr 


2 rs Shae hee been. assumed to ‘behave as real gases. . It ‘is of’ interest to see ene 


shes titeg  Peketete NE ha. “vhat would. happen to the therinodynanie properties of the propellant,” co oe 
ae 2 ee Pe ee sci Sa 
BAD beg eee a oe if this assuniption were: not martes eg., assume: the combustion products 


behave as ideal gases. In “the; ‘Brinkley code this ‘is done by setting. : 


a ae pk t _ ‘the virial coefficients equal to er. In the appendix: are listed. - 
4. ets . 


te my wee Ut the thernodynante ‘properties for ML propellant. at a density of load~ 
R eS Sha as _ ing of 0.2 gm/ec. Comparison of t hermodynamic propenties. ae one. 
oe bags EA 9 wT pe ae ‘temperature ervey the following results: 


a cS ae ae . aay ‘e 23 a . mt “ “yh a 


Nba Arties Prd ete aI Cele Ae 
secant Aisa of STDS DAT a ahem I Reals ha 


- Real Gos, ae ine “Taeal Gas :. see Pe oe 2) 
“ A B35 a re - aoe ee ee 


~ dn the ‘constant. ‘voluine ectitastion. préseure, the facet gas. cénbustion oo aoe 
| Peas. being: 2h 38 lower ‘than’ the real, gas conbustion: bressure.- 


_ THIERYODINANTG DATA 


: hermodyianic data ‘used. tn “ene: ‘eouputetio ons are. préaented fa’. 


, “ables ok through’ ‘Tab?e Xaite ‘Table. sag ‘shows. the atomte composition: . 


ae and, heat of formation of “wil itary: and exper ize tal propellante, Thesé. 


ie “data, ‘cogether vith ‘the. iadting: density, are: all that ‘are. needed for. 


e ‘propellant “input to ‘the Brinkicy program. ‘Tr. ‘Pable x “the atomic’: oo oe 
ae composition and. heats oF. formations’ are. given for “tlie ‘special’ pron 


_ pellants, ‘that ‘is, 2 “the ma propellant containing various’ “eoncentra~ ms 
tions of. ‘water. and ethyl slcohol, and th. experiiiental dovble ‘pases 
as “propellants Table XI nhows the atomia compositici.. and. heats: of 


“fovmation for the propellant ‘reactants. Yucluded in’ chi. table ie 


se ‘the souitce of the ‘thermedyiante fate a used to. compue the heats. of 
Coes a . m - formation of the’ Feact ants. - ca eae eae : wee 
x e - co : ee ie e Ns it ‘was ‘absumed ‘iy compitting the heats ‘of t ‘Foiu-bion oft the ‘pio= ~ a 
: 2 . aay oe : = e on : pellente considered: in ‘this report that ‘the ‘total: ‘hext of formation” . 
es te ; o ar “Was, the ‘sum of the proportional heats of. formations of phe individual 

a oa i. os “ - reactants. Thus heats of polymer}. zation, ur auy heat released by’ : 
Bee eR. ENT a ‘chemical reaction between the propellant reactants prior’ te e burning, - 

sae not. coneidsred. x 


a ti guetta eRe emis ‘2 “the teats: of Pormatibe: of the: propellant cambust t ton. products are - 
. Pe gles ee given’ cCy ‘Table XIL F wand in Table XIII, a correction factor to correct" 
ft 4 Oe . : “ ; : , : 
NG rey 

ty 

a 

a 15 


ae] 
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, ULES AT ES 


fei “the vais of ‘the. taternad ‘energy ‘of. the combustion peices from rs 
ena : Level, ‘of o%k (the’ level being ‘used in. the Brinkley caldulations) to 


“8 ‘evel of. 298° Ky the temperature’ level or the: heats of formaton. . 
a of: the: propellant reactents, 


‘Tables. of internal energies, cheintoal.: 
‘potentiat, “specific heats, and virial, curtection ‘for edch of tie ee 


combusts on. Products are presented ' ‘fini Lesert s 2 reyort(2) 
: be Fepeated heres . < : ae 


ee RAGE va SORT 


and will, not 


hs ¥ DDISCUSSEON 


oe thie’ ‘report presents thiermodyiicinte' data, "veaeiding* the ‘combustion | 
"properties of military gun ‘propellants, plus. some data concerning . 


-expertnisatel, propellante. ; Phase data: must be regatded rCE “tentative, - 


_ pending: revisions i the. thermodynamic date, used in the computations. 
ah “In future’ vevi sions," doth | the, ideal gas therniodyniante properties, 


LEO pg 2 


ee 


and. the. ‘real: Be. S$ vi rial ‘corrections of ‘the combustion products will 
“be: replaced. ‘In ‘adgstion, the thermodynante: code will ve set “up in.’ 


; “sudh a way! as to ‘account tor. the yresencé. of the tnorgante sata’ 
“aise. dn the propeLlants.: os we 


Ta Leser's seport (3) I, ‘one serious error occurred ‘unten, until 


x “ Loéated ‘and cor “ected, markedly “influenced the shape OPthe oo oS 
_pressure-entropy. diagram,’ 


s 


This’ error was’ ia the ‘sign. of the. 
‘correction’ ‘to. the’ activity coefftetent for gas tipartection. 


the . 
” eorrest ‘equat ton is [en (ay an Leser's hae 


. #306 7, « gcEee ele. soe 1/2 et 


a ae oe ea a , ‘the ‘ideal gas therinddynianté data and the real gas virial | 
Bat es ge data for amon had aiscontinuties and drregularities when ° 


ee plotted as @ ftinction ‘of ‘temperature, ° In Fig. 15, the effects of 


ee _ ‘these irregularities apon the combustion gas composition for MeL ; 
ME at 0,2 gn/se donatne. is shown, Since the proportion of ammonia in 


‘ “Ps ak ote fe ; Se . . x 16 


et e oY 
we tt 
x ta) 
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ni - are gay —arerr ss _ oy > ae = — = ae ae 9 om pose nals 
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ae ies eonbustion 2 products. is’ very minor) ‘any changes ‘tin the ‘asinonte 

: data, would hot be: ‘expecticd ‘to change. the overall thermodynamic 

a ve properties. of the propellant - to any great ‘extent. A severe error . 
Sg in the chemical potential. for eumonia at 1700°K wes ‘corrected in 

es “the. process: of these commutations: an is not (Present an the resulte, 
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5 OF TABLES | 


Ba anes < 7 ae Hie : “Gn the following tables thernodynamic properties end combustion 
an . ' a : gas. ‘composition are listed for all the propellants which appear an 
Tables I through IV. In the thermodynamic properties and combustion 
ges: composition tables): the type or propellant, for instance Ml, sppears - 
“4m the upper left of the title block. “In the upper right of the title 
oS is babe the loading density ; (gas density); £ for instance, 0.05 gnfec. 


“Tn the row feLlowing temp. °K ace combustion gas tenperatures at 
which the combustion g22 thermodynamic properties and composition wes 
"computed. These temperatures are ‘Usted starting ata temperature 
considerably below the flare temperature, and inerease in 100 eK. 
“| intervals until a ‘temperature is reached which is usually about 
“sre 200% K above ‘the flame temperature. The last ‘temperature in this 


cg ee taRE ie ‘the ahaes venperatures 


“To the column dlow each temerature are the combustion gas 
| thermodynamic properties and composition at a temperature and gas 
~ density (Losiding beater The thermodynamic properties are ss 


_ follows: , 


- NOL.WT. - Average woleculer weight - gmn/gm - mole 


RS Pe _ = Pressure - psi 
He e _~ Energy content - cal/gm 

Bee @. 
eo cv -.= Specific heat at constant voluse - cal/gm - “K 
a4 y - Ratio of specific heats , 
i : , - Oi: 

ee S - Entropy - cal/gm - K 

e - Internal energy - cal/gm 
HF - Heat of formetion - cal/g. 

i! The thermodynamic values are complete for all temperatures 
See : i ‘ 

aa listed, exeept the flame cemperatiure., In the Clame temperature 

Me Oy . . 
a . oe eoluwon the velue for energy content was not compuled, thus a black 
Bo . 

x : ES os space appears for this value, 
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_ Below the tublé of thermudynanic proverttes are listed the 
combustion gas ccmpesition. The gas spectes are: CO, C05. Has 
#0, Ny» 0,, 0, OH, H, NO, H, iH» and CH). Th2 units of compoesi- 
tion are gn - moles/gn. — 


In the gas composition tables e number of vlanks appear. - The 


presence of blanks means that the computed concentzatioa of a parti- : 


Bh : : o g cular specie was less then 1 x 107! gm - moles/em for the tempera- 
aaa ; , ; 


ture and gas density shown. 
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-- CONPURATION OF PROPELLAN? FLAME TEMPERATURE 
(3 y. 


Tr Leser' s veport ne shouiston was made in ‘the code a the. aoe 


* computation of the propellant's flame. ‘temperat ture. The. following 
a: describes the method used in the revised code ae) that the flame 
temperature of the propellant could be computed. sy eae: 


Tn sant) a ‘method ic given by which the flane temperature of: a ane 
- propellant may te computed, In this reference the ‘flame temperotur2.” 


ae aefined as the temperature’ at which the ‘internal euergy of the 


combustion gases is equal. ‘to the heat: of reaction. oF the propellant. re ae 


This can ‘te “epresented asi e 


where: | | 


AHF = HF prod - WF prop. = heat. of reuctin of propellant. - 


en on oe internal energy of conbustioa | gas ‘at flame temperature. = 


cal./em 
HF prod = heat ae of ‘combustion gases - beaten 
_-HF prep = heat of formation of propellact . “ . cal/em 
The heat of formation of the Propel lant is: given for each pro- 
pellent, and as such. zesd dato che program as a constant. : The 
‘dase temperature used to compute heats of formation is 25 ie (298°K). 


. For come propellants heaty oi formation were given at 300%K.. The 


difference in heate of formations between the two (emperatures was 
assumed to be negiigible. 


The heat of formation of the combustion gasea is a function cof 
the combustion gas composition. so il has to be computed at each 


trial temperature for which comtustion gas composition is known, | 


“The function used ie: 


a 


WF prod = & 9, (AH, £ (2 =1,2,-~- 8) (2a) 
iel 


cai! ee number at moles of isth ‘eonbustiion product - + nees/i 


(an, de: 


Rent, of: formetion of i-th combustion produst-tal/s ole: 


ie. 


: ‘the. heats, of formation. of the con ways tion 2 products are given 3 ‘in : oe 
tents XID, + ae: : 


fad. ‘the: orig sinell Brinkley cote, an” energy conten’ "e 1 is” computed: 


en is defined ‘as. tne energy. of, the. coxbustion gases at erature ee 
a, relative ‘to. ‘the ‘elenents at 0 Ke Sthea the: internal energy ‘of the, 
7 gases Pigs is defined es the “ainount, of endrgy required to raise the » 
" combustion gas 8 from SOE datum vébperatore: yin this case 29. °c) 
ke some temperature T; it is ‘spparent vhetl a. correction, ‘factor met 
“be used to convert ‘the. compute energy ceatent * a" lo tte faternal » 


. energy mM “ths conversion As. nade by ‘means of the following: 


relation: °°; ; beet ae _ on mee an 
ee ee eer : pa Sate Be 
ee = ee > n, 8 # “2 hoes She an 
+ . £ L oad 1° : - " : ie ie 
“where: Bog _ “change ha internal, ene erEY ‘of ieth  conbustion oe 


‘product from 0° K to est Ke: .cal/nole 
heat of ‘formation of i-th conbustton peoduet 


a) 


i ‘from the eleneits: at OK - eel /mote | 
‘Tables of, these, vals zes for all the coubustion ‘products ‘is. found ~ 
in NBS Circular 500." n Table xIIT, ‘the correction, factor a 


As listed for ea: n of ‘the conbus stion products. : 


axe) determine ‘the ‘Plane temperature, the computer prograun 
computes E, HF prod, AUF, end APF- for a temperature range. which” 
has been previously chosen so as to be several huadred of degrees 
“above and below the estimated flere temcerature = upon completion 
_of this portion of the problen, the conputer checks beck until it 
reaches a set of temperetures in wittel vhe difference AHF-b has 
changed from a positive value to a negative value. Using a 3 
point Lagrangian Internolation forculas on values o7 the 
ditYereuce at three te nparature Steps; “790 of which ere positive 


and one negative, a temperature value {£3 found in wich thts 


at 


Septet ee Minti = tne ert me oe 


EERE SE rete ae SRE ENTE Fee ae AN SRNR REE oe RIS fahemmmeaeS aptamer 


Peta kn BE ples 5G Got Nae 


difference ‘Ls oe ‘This tex perature: is ‘the ‘lone temperature. - eee 

In addition,. the inte rpolatiod. ratio are used ‘to. “eompute ‘thermo+ ae 

a Aynamic. and gas compos tion. val 128 at the fl eine temperature... When = 
“this” port ion of the program ts. conpleted, : 


‘the: values ure 2 printed 
2. cout ‘end ‘the ) Progran: is. atopped. ere ase " ‘ nit 
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